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SUMMARY 

•When  the  convection  velocity  of  the  turbulence  exceeds  the 
local  speed  of  sound  in  the  medium  the  coupling  of  the  panel  modes 
to  the  acoustic  field  becomes  very  strong.  The  structural  damping 
which  is  usually  constant  for  all  modes  is  overshadowed  by  the 
acoustic  damping  which  increases  with  frequency.  For  a  boundary 
layer  of  nearly  zero  pressure  gradient  the  linear  theory  is  still 
useful  in  predicting  the  panel  response  and  compares  well  with 
experiment. 

Using  Kirchhoff's  theorem  the  acoustic  power  radiated  from 
a  turbulent  boundary  layer  was  expressed  in  terms  of  the  cross 
spectral  density  of  the  pressure  field  over  a  rigid  plane  surface. 
Computation  shows  that  by  using  this  method  the  acoustic  power 
level  agrees  satisfactorily  with  the  measurement. 


I.  Introduction 

A  large  aircraft  in  supersonic  flight  opens  new  challenges  in 
aeroacoustlcs,  the  pursuit  of  which  is  vital  to  advancement  in  this 
field.  Problems  emerge  because  subsonic  experiences  are  Inadequate 
and  it  becomes  increasingly  difficult  to  define  a  clear-cut  analysis 
for  predicting  the  response  of  a  panel  to  a  flow  field  on  an  aircraft 
fuselage.  The  aerothermoacoustic  and  aerothermoelasticity  problems 
are  not  part  of  this  investigation.  In  the  present  experiment  the 
wall  temperature  averages  around  550  degrees  Rankine  and  the  response 
of  the  panel  is  linear,  no  stress  deterioration  or  stiffness  change 
will  occur.  Of  course,  this  is  a  specialized  case  of  nearly  zero 
pressure  gradient  where  the  power  spectral  density  technique  of 
generalized  harmonic  analysis  can  be  used.  In  the  case  of  larger 
pressure  gradients  (shock  wave)  a  more  rigorous  approach  is  required, 
since  the  flow  field  leads  to  Instability  (flutter)  Dowell  (1966-1967), 
Maestrello  (1968).  It  then  becomes  necessary  to  use  nonlinear  plate 
theory  and  convenient  to  use  time  rather  than  frequency  domain  analysis. 
Even  with  the  simplifications  of  linear  analysis  and  use  of  the  power 
spectral  density  technique,  the  problem  is  still  complicated  due  to 
the  necessary  inclusion  of  the  pressure  resulting  from  panel  motion. 

The  forcing  field  terms  excluding  the  pressurization  effect  include 
the  pressure  over  a  rigid  wall  and  that  produced  by  panel  motion. 

The  resulting  quadratic  and  inhomogeneity  gives  rise  to  a  reflection 
coefficient  which  imposes  additional  boundary  conditions. 
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This  last  development  is  contained  in  a  report  which  will  follow  by 
Maestreilo  and  McDaid.  The  present  report  is  concerned  with  obtaining 
radiation  from  a  turbulen  boundary  and  panel  response  to  it  with 
emphasis  on  the  coupling  of  panel  modes  to  the  acoustic  field  which  is 
predominant  when  the  turbulent  convection  velocity  exceeds  the  local 
speed  of  sound  in  the  medium. 

These  measurements  were  made  on  the  side  wall  of  the  Jet  Propulsion 
Laboratory  20  inch  supersonic  wind  tunnel.  The  boundary  layer  is  turbulent 
on  the  side  wall  from  natural  transition.  The  tunnel  was  well  suited  for 
these  measurements  partially  because  the  flow  conditions  are  well 
documented  (Coles  1963)  and  partially  because  the  statistical  measurements 
of  the  pressure  and  velocity  field  are  well  recognized  (Laufer  1961, 

Kistler  and  Chen  1962). 
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II.  Experimental  Arrangement  and  Instrumentation 

The  flow  investigated  was  the  side  wall  boundary  layer 
of  the  JPL  20  inch  supersonic  wind  tunnel.  The  tunnel  is  a  continuous 
flow  type  utilizing  a  two  dimensional  flexible  nozzle  which  provides 
for  speeds  up  to  Mach  number  5.  A  description  of  the  facility  with 
a  range  of  operating  conditions  is  discussed  by  Laumann  (1967). 

For  a  brief  description  with  r, articular  emphasis  on  those  features 
which  are  most  important  for  velocity  and  pressure  field  measurements 
one  may  consult  Laufer  (1961)  and  Kistler-Chen  (1962) . 

The  experiment  was  arranged  to  perform  three  basic  measurements 
mean  velocity  profile,  wall  pressure  fluctuations  and  the  response  of 
a  simple  panel  structure.  The  side  wall  window  was  modified  to 
accommodate  two  identical  rigid  steel  plates  which  supported  the 
instrumentation  required  for  measurements.  One  plate  contained  an 
array  of  holes  to  mount  the  pressure  transducers.  Figure  1.  Another 
had  a  cut-out  portion  in  which  the  test  panel  was  mounted.  Figure  2. 

The  panel  was  tested  with  a  static  pressure  differential 
loading  of  0.06  and  14.0  psi.  The  operating  static  pressure  of  the 
test  section  is  always  below  atmospheric  pressure.  In  the  case  of 
0.06  psi  pressure  differential,  a  chamber  of  size  8  x  14  x  6.6  inches 
was  built  around  the  panel  to  contain  the  pressure  which  was  kept 
equal  by  interconnecting  the  chamber  with  the  tunnel  by  a  row  of  holes 
located  downstream  from  the  panel,  Figures  (2)  and  (3).  Tests  with  the 
14  psi  differential  were  run  both  with  and  without  this  cover. 


|  Displacement  Probe 


rr 


1 


( 

* 

i 

i 

\ 

i 


! 

I 


,«>" 

-* 


The  boundary  layer  on  the  tunnel  side  wall  is  turbulent 
from  a  natural  transition  on  a  smooth  surface.  Measurements  were 
made  in  the  region  of  zero  pressure  gradient,  as  well  as  under 
conditions  of  shock  and  boundary  layer  interaction.  The  shock  was 
induced  by  a  30°  wedge  mounted  outside  the  boundary  layer,  off  center 
and  to  the  same  sice  where  the  measurements  were  made.  This  was  done 
to  offset  the  position  of  the  reflected  shock  from  the  opposite  side. 

The  position  of  the  shock  was  determined  by  observing  the  displacement 
of  a  line  of  tufts,  and  more  precisely,  by  a  static  pressure  survey. 

The  response  of  the  shock  interaction  with  the  panel  will  be  reported 
separately. 

The  velocity  profile  was  measured  by  a  standard  technique. 

The  wall  pressure  fluctuations  were  measured  by  pressure  transducers 
mounted  in  the  plate  shown  in  Figure  1.  The  transducers  were  mounted 
on  a  threaded  bracket  in  order  to  adjust  the  alignment  within  ±.0005  inch. 
The  alignment  was  checked  with  a  micrometer  at  each  transducer  location 
at  the  tunnel  operating  pressure.  This  was  the  most  time  consuming 
operation  of  the  test.  If  the  flushness  is  not  maintained  within  the 
tolerance  the  power  spectra  would  assume  a  different  shape  in  the 
lower  frequencies. 

Two  types  of  pressure  transducers  were  used;  one,  the 
conventional  lead  zirconate  titanate  made  by  Atlantic  Research  and 
the  other  a  capacitance  type  made  by  Photocon  Corporation  with  a 
sensitive  diameter  of  .060  and  .090  inches  respectively.  This  last 
type  of  transducer  was  used  to  survey  the  low  frequency  part  of 
the  spectrum  because  it  is  less  sensitive  to  vibration  than  the 
crystal  type.  A  monitoring  survey  of  the  power  spectra  and  cross 
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correlation  was  made  during  the  test  to  insure  that  satisfactory 
measurements  were  recorded  on  the  magnetic  tape.  In  line  power 
spectral  density  and  cross  correlation  were  obtained  with  a  Hewlett 
Packard  wave  analyzer  and  with  a  Boeing  built  correlator.  The  output 
was  recorded  on  an  Aapex  FR-180QH  14  channel  tape  recorder  in  the 
FM  mode.  The  main  feature  of  this  recorder  is  that  it  has  a  low 
static  and  dynamic  skewness.  Four  channels  were  used  for  simultaneously 
recording  data  for  correlation  measurements,  two  channels  for  the  lower 
frequencies  down  to  the  dynamic  range  of  the  recorder,  and  two  channels 
for  the  higher  frequencies.  The  maximum  dynamic  range  was  obtained  by 
splitting  each  data  channel  into  two  tape  tracks  through  phase  matched 
filters  to  separate  the  lower  and  higher  frequencies.  A  high  frequency 
sine  wave  reference  tone  was  recorded  on  alternating  tape  tracks 
simultaneously  with  the  data  for  the  purpose  of  correcting  time  error 
between  data  channels  used  for  digital  correlation  analysis.  This 
technique  provides  a  broad  band  correlation  in  both  frequency  and  dynamic 
range. 

The  test  panel  is  thin,  rectangular,  titanium*  of  size 
12  x  6  x  .062  inches  brazed  on  all  four  sides  to  a  3/4"  x  3/4" 
titanium  frame.  The  brazing  is  intended  to  simulate  the  clamped 
edge  condition.  The  frame  and  skin  assembly  was  placed  in  a  vacuum 
retort  using  Handy  and  Harmon  braze  alloy**. 


*  T1-6AL-4V  Titanium  alloy  containing  6%  Aluminum,  4%  Vanadium 
and  90%  titanium. 

**Braze  alloy  70%  Ag,  28%  Cu,  2%  Li. 


The  retort  was  placed  in  a  furnace  and  raised  to  a  temperature  of 
1450°F.  The  panel  assembly  was  bolted,  into  the  cut-out  portion  of 
the  plate  and  made  flush  with  it.  The  panel  displacement  was 
measured  with  "Photocon”  capacitance  type  displacement  transducers 
mounted  on  brackets  which  could  slide  along  a  bar  and  could  be  set 
precisely  by  means  of  a  screw  mechanism.  The  transducer  mounting 
can  rotate  360°  to  allow  correlation  along  the  longitudinal  and 
lateral  axis.  In  line  power  spectra  and  cross-correlation  were 
taken  and  then  recorded  on  magnetic  tape  for  0.06  and  14  psi  pressure 
differentials.  Measurements  were  repeated  with  damping  tape  placed 
on  the  outside  surface  of  the  panel.* 


o 

*Two  layers  of  damping  tape  with  surface  density  of  one  layer  0.112  Ib/FT  . 
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III.  Measurements  and  Discussion 

Three  basic  measurements  were  made:  mean  velocity  profile,  wall 
pressure  fluctuation,  and  the  displacement  response  of  a  simple  titanium 
panel  with  clamped  edges.  Wind  tunnels  inherently  have  a  certain  amount 
of  interference  (vibration  and  sound  propagated  up  or  down  stream 
through  the  boundary  layer).  Commonly  the  vibration  interference  can 
be  bypassed  on  the  measurement  of  the  wall  pressure  using  a  capacitance- 
type  transducer.  In  the  present  test  measurements  of  the  wall  pressure 
fluctuation  were  rejected  below  800  hz  due  to  acoustic  interference. 

The  tunnel  vibration  will  not  alter  the  dy“:jnics  of  the  panel, 
since  it  is  small  compared  with  that  induced  by  the  boundary  layer. 

The  fundamental  mode  (500  and  600  hz  for  the  unpressurized  and  pressurized 
case),  will  be  effected  by  the  tunnel  interference.  The  effect  is 
predominant  on  the  fundamental  mode  and  is  not  expected  to  be  larger 
than  a  factor  of  two. 

a)  Mean  Velocity  Profile 

The  quantities  directly  measured  were  the  impact  pressure, 
wall  pressure  and  stagnation  temperature.  The  stagnation  temperature 
and  static  pressure  are  assumed  to  be  constant  throughout  the 
boundary  layer.  The  local  static  temperature,  density  and  velocities 
are  obtained  from  the  energy  and  state  equations.  Therefore,  one 
can  assume  that  the  properties  of  the  boundary  layer  closely 
approximate  the  properties  of  the  equilibrium  adiabatic  flat  plate 
boundary  layer. 

Measurements  made  by  Coles  (1963)  using  the  above  assumption 
over  a  flat  plate,  estimated  an  error  in  momentum  thickness  6 
of  about  2%  up  to  Mach  5.  Table  1  lists  some  of  the  most  pertinent 
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flow  parameters  in  a  supersonic  flow  at  zero  pressure  gradient.  The 
mean  velocity  profiles  at  various  Mach  numbers  were  reduced  to  a 
constant  density  form.  The  method  used  is  based  on  the  work  of 
Coles  (1963)  in  which  he  assumed  that  in  the  sublayer  the 
Reynolds  number  remained  unaffected  by  the  compressibility  provided 
that  the  density  and  viscosity  are  evaluated  at  mean  sublayer 
temperature.  Ke  proposed  a  simple  generalization  of  the  law  of 
the  wall  such  that  the  similarity  functions  remained  the  same  as 
in  incompressible  flow.  The  transformed  law  of  the  wall  is  obtained 
by  adopting  the  Dorodnitsyn  type  scaling  which  produces  a  lateral 
stretch  of  the  coordinate  providing  an  equivalent  incompressible 
profile.  This  transformation  coincides  with  the  Von  Karman 
hypothesis  that  the  skin  friction  can  be  obtained  from  the  in¬ 
compressible  relations  by  evaluating  the  density  and  viscosity  of 
the  wall  instead  of  the  free  stream.  An  extension  of  Coles  approach 
has  been  recently  made  by  Baronti  and  Libby  (1966)  and  used  by 
Watson  and  Cary  (1967). 

Their  conclusion  was  that  for  an  adiabatic  flow  up  to  Mach  5 
the  Dorodnitsyn  type  transformation  can  be  used  to  represent  the 
boundary  condition  such  that  a  universal  description  of  the  law  of  the 
wall  is  obtained,  but  above  Mach  5  their  method  was  not  successful. 


From  the  law  of  the  wall  in  the  logarithmic  region 
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U_  =  A  In  B  5  ?f 


(1) 


where 


dy. 


(2) 


and  in  the  laminar  sublayer 


U_ 

U 


0  <  C  i  L 


(3) 


i,  _  ry 

0/0  *  U/U  *  (C,/2)^f(C)  and  I  p/p  dy  is  the  Dorodnitsyn 

T  e  f  JQ  e 

_  _  u 

variable,  U  *  (t  /p)  ,  C£  *  skin  friction  coefficient,  y  =«=  coefficient 
t  oj  t 

of  viscosity  A  =  2.43,  B  =  7.5  and  *  10.6.  The  ratio  yec/u  comes 
from  the  assumption  due  to  Donaldson  (1952)  that  the  Reynolds  number 
associated  with  the  laminar  sublayer  is  invariant.  The  edge  of  the 
boundary  layer  was  selected  at  a  point  where  the  ratio  U/Ue  *  .995.  The 
ratio  p/pe  is  obtained  by  using  Crocco's  relationship  for  the  density 
profile.  The  compressible  skin  friction  coefficient  is  related  to 

the  incompressible  by  the  following: 


^f  =  (peUe/pe,lia,)  (w/owe)Cf 


(A) 


The  values  of  Cc,  Cr,  y  o/y  and  x  are  given  in  Table  1. 
f  f  e  a) 

A  plot  of  the  law  of  the  wall,  Figure  A,  using  the  above 


WWW6G3 


15 


transformation  compares  well  with  the  measurements  of  Laufer 
and  Klebanoff  which  were  made  for  an  incompressible  flow.  It  can 
be  seen  further  that  the  present  data  transform  for  an  equivalent 
incompressible  boundary  layer  is  in  close  agreement  with  the 
relationship  found  by  Coles,  Figure  (5),  as  well  as  with  measure¬ 
ments  made  by  other  investigators. 


Fig.  5.  Local  skin  friction  coefficient  for  equivalent  incompressible 


b)  Wall  Pressure  Fluctuation 


The  problem  is  to  obtain  an  appropriate  description  of  the  wall 
pressure  fluctuation  produced  by  a  turbulent  boundary  layer  at  near  zero 
pressure  gradient.  From  a  spatial  and  temporal  correlation  of  the 
pressure  field  an  estimate  of  the  response  of  the  structure  can  be 
made. 

The  structure  senses  the  different  components  of  the  pressure 
as  they  are  swept  past  the  surface  at  different  velocities.  Since  the 
panel  response  behaves  like  a  narrow  band  filter,  a  meaningful  represen¬ 
tation  of  the  field  can  be  obtained  with  less  complexity  than  would  be 
required  to  represent  the  turbulent  flow  itself.  The  convected  nature  of 
the  wall  pressure  fluctuation  can  be  described  by  a  space-time  correla¬ 
tion  function.  The  convected  feature  is  not  independent  of  either 
spatial  or  temporal  delay  on  which  they  are  measured  but  is  contingent 
on  both.  The  decoupled  equations  for  a  compressible  viscous,  heat 
conducting  gas  within  the  framework  of  the  linear  theory  have  three  types 
of  disturbance  fields.  The  predominent  ones  are  the  vorticity  and  the 
entropy  modes  whose  phase  velocity  is  the  usual  convection  velocity.  The 
third  type  is  the  acoustic  mode  in  which  energy  is  transported  away  at 
the  acoustic  wave  speed.  All  three  modes  are  independent  except  at  the 
wall  because  they  are  imposed  on  the  velocity  and  temperature  (Kovasznay 
and  Chu  1958). 

Laufer  (1964)  successfully  measured  the  acoustic  contribution  by 
recording  the  density  fluctuations  along  the  Mach  wave  outside  the  turbulent 
boundary  layer  using  a  hot  wire  anemometer.  Wilson  (1959)  measured  the 
near  and  far  field  of  noise  radiated  by  a  turbulent  boundary  layer  developed 
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on  a  rotating  cylinder.  Separating  the  contribution  of  the  acoustic 

mode  at  the  wall  from  the  other  disturbances  is  more  cumbersome  because 

the  energy  radiated  away,  from  the  turbulence  region  represents  only  a 

small  fraction  of  the  kinetic  energy  flux  of  the  turbulence  and  an  even 

smaller  fraction  when  compared  with  viscous  dissipation.  However,  the 

wall  pressure  fluctuation  measurement  retained  the  contribution  of  the  acoustic 

mode  as  well,  and  the  non-zero  correlation  integral  can  be  interpreted  as 

energy  dissipated  acoustically.  The  sound  sources  are  located  away  from 

the  surface  of  higher  order  than  dipoles  whose  strength  is  zero  at  the 

surface.  The  sources  are  equivalent  to  a  quadrupole  system.  The  surface, 

assumed  to  be  rigid  and  large,  acts  as  a  simple  reflector  of  sound, 

Doak  (1960),  Ffowcs  Williams  and  Lyon  (1963),  Lyams’nev  (1961),  Powell  (1960), 

Meecham  (1968),  and  Kraichnan  (1957),  apart  from  the  fact  that  a  contribution 
to  the  pressure  may  occur  due  to  the  viscous  effect,  Leehey  (1968). 

Measurements  of  the  wall  pressure  fluctuation  were  made  in  the 
range  between  0.8  to  200  KH^  using  the  conventional  lead  zirconate  ceramic 
disk  as  well  as  with  capacitance  transducers  with  diameters  of  0.060 
and  0.090  inch  respectively.  The  normalized  power  spectra,  corrected  for 
attenuation  due  to  the  finite  size  of  the  transducer  (Corcos  correction),  are 
shown  in  Figure  6.  The  deviation  in  the  shape  at  lower  frequencies  may  be 
attributed  to  extraneous  interference*  which  persisted  inside  the  tunnel  and  to 
radiation  from  the  turbulence  (this  last  one  is  the  apparent  surface  source 
originated  by  the  source  near  the  surface).  The  extraneous  interference, 

Narrow  band  spatial  correlation  measurement  for  frequency  below  0.8  KH^ 
using  capacitance  probe  transducer  (quite  insensitive  to  vibration) 
shows  that  the  decay  of  the  turbulent  pressure  has  superimposed  an 
acoustic  component  similar  to  that  shown  by  Wills  (1968). 
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Fig.  10.  Space  time  correlation  of  the  wall  pressure  fluctuation. 
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Fig.  11.  Space  time  correlation  of  the  wall  pressure  fluctuation 
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however,  masks  the  radiated  noise  from  the  turbulence;  this  is  the  rea¬ 
son  that  the  frequency  spectra  have  been  truncated  at  0.8  KH  .  Lilley 

z 

and  Hodgson  (1960)  predict  that  the  wall  pressure  spectral  density  varies 
like  the  square  of  the  frequency  at  low  frequencies  for  vanishingly 
small  Mach  numbers.  Ffowcs  Williams  (1965)  shows  that  for  finite  com¬ 
pressibility  the  correlation  area  varies  like  the  square  of  the  Mach 
number.  The  present  model  of  the  wall  pressure  fluctuation  retains  some 
of  the  features  contained  in  the  above  comments . 

The  normalized  power  spectral  density  (Fig.  6)  is  similar  to 
that  obtained  from  subsonic  flow  over  a  flat  plate.  It  peaks  at 
t»6/Ue  ~  2.0,  Bull  (1967),  Serafin  (1963),  and  Hodgson  (1962),  show 
similar  behavior  (Fig.  7).  In  contrast  for  pipe  and  duct  flow  the  wall 
pressure  spectral  density  indicates  no  predominant  peak  Cor cos  (1960), 
Maestrello  (1965).  In  flight,  measurements  made  of  the  wall  pressure 
fluctuation  have  a  more  difficult  interpretation  because  in  general  there  is 
an  apparent  pressure  gradient,  as  well  as  some  roughness.  Hodgson's  measure¬ 
ment,  however,  over  the  glider  wing  produced  the  typical  rolling  off  spectrum 
at  lower  frequency  while  measurement  made  on  the  727  airplane  at  higher 
subsonic  Mach  number  indicates  no  significant  rolling  off  at  the  lowest 
frequencies,  Maestrello  (1967).  Schloemer  (1966)  from  a  set  of  systematic 
measurements  points  out  that  the  lower  frequency  part  of  the  spectral  density 
is  sensitive  to  the  pressure  gradient.  A  more  extreme  case,  that  is  for 
separated  flow  Maestrello  (1968)  indicates  a  similar  observation. 

Figs. (6  and  7)show  two  main  eddy  structures,  the  larger  scale,  faster 
moving,  lower  frequency  size  structure  scaled  with  the  boundary  layer  thickness, 
and  the  smaller  scale  higher  frequency, slow  moving  eddy,  scaled  with  the  sub¬ 
layer  thickness.  This  is  also  born  out  from  the  law  of  the  wall  by  application 
of  the  transform  for  constant  density  flow  discussed  in  the  previous  section. 
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Further  results  show  that  the  ratio  of  the  wall  pressure  with  the 

local  wall  shear  stress  p'/t  is  nearly  independent  of  Mach  number,  while 

w 

the  ratio  of  the  wall  pressure  with  the  free  stream  dynamic  pressure 

p'/qe  ~  R^.  The  results  are  reported  in  Figs.  (8  and  9).  Similar  results  have 

been  previously  obtained  by  Kistler  and  Chen  from  measurements  made  in  the 

same  wind  tunnel  except  for  about  30%  discrepancy  in  the  values  of  p’.  From 

the  figures  the  product  p'/qe  *  t  /p*  reproduces  the  well-known  results 

of  the  frictional  law,  the  compressible  skin  friction  coefficient  decreases 

as  the  Mach  number  increases  at  a  constant  Reynolds  number. 

The  broad  band  space-time  correlation  curve  for  separation  in 

the  stream  direction  R(?,0;t)  and  in  the  lateral  direction  R(0,n;t)  are 

shown  in  Figs.  (10  and  11),  obtained  at  constant  spatial  separation  and  at 

variable  time  delay.  The  cross  correlation  retained  a  frequency  band-width 

from  0.5  KH^  to  in  excess  of  100  KHz  with  an  adequate  dynamic  range. 

The  longitudinal  correlation  shows  the  usual  convected  feature,  slowly 

losing  its  coherence,  while  with  the  lateral  correlation  no  convection  is 

observed  except  possibly  at  large  separations  where  the  envelope  flares  out 

at  small  correlation  values.  The  broad  band  convection  velocity  corresponds 

to  U  /U  =  0.78  at  M  =  1.98  and  U  /U  =  0.72  at  M  =  3.03.  It  is 
c  e  e  c  e  e 

interesting  to  point  out  that  for  correlation  measurements  taken  at  constant 
spatial  separation  at  a  slight  angle  off  the  flow  direction  the  convection 
velocity  was  slightly  greater  than  the  one  obtained  along  the  flow  direction. 
This  could  possibly  be  explained  as  disturbances  moving  with  particular 
orientation  for  the  case  of  supersonic  convection,  since  it  is  a  propagating 
wave  system.  Unfortunately  insufficient  data  was  obtained  to  construct 
the  wave  pattern  over  the  surface. 


26 


•  ?T?s*W.jfr  .  '■■■  '  -v  - 


From  measurements  made  of  the  wall  pressure  fluctuations,  plotted 
dimensionally  in  Fig.  (6)  a  model  is  obtained.  The  model  is  of  similar  form 
to  that  reported  by  this  author  earlier  (1967  and  1968)  for  a  duct  flow 
except  it  takes  into  consideration  a  peaked  power  spectrum,  characteristic 
of  a  boundary  layer  developed  over  a  flat  smooth  plate  at  zero  pressure 
gradient.  From  the  figure  the  power  spectrum  P(u)  and  the  normalized 
spectrum  7i(m)  satisfy  the  following  relationships: 


dm*  1  J  P(w)dm  *  <P2>  ss  12t2  ; 


Tt(a))U  4  -K  (oj6/U  )  (5) 

_ £  =  Y  A  e  n  e 

6  4  n 

n=l 


where  *  4.4  x  lO-' 
A2  *  7.5  x  10  ; 

*  -9.3  x  10 
A,  =  -2.5  x  10' 

*4 


Kx  =  .5.78  x  10 
K2  *  2.43  x  10' 
K3  =  1.12 
K4  =  11.57 


The  normalized  auto  correlation  is 


r(t>  =  y 


,  K2  +  (tU  /S) 2 
n=l  n  e 


The  flow  considered  is  semi-frozen,  that  is  with  decay  in  space  and  time 
at  constant  convection  velocity  Uc.  From  Eq.  (6)  using  Taylor's  hypo¬ 
thesis  the  longitudinal  spatial  correlation  is: 


R(0  -  J 


~  /5U  \2 

n-l  K2  +  hfr 

n  \  ou  / 

\  c/ 


(7) 


From  experimental  observation  the  lateral  spatial  correlation  suggests 


R(n)  =  e 


■n/a2(S 


and  the  normalized  cross  correlation 

4  -Ui/0^  -|n|/o26 

A  K  e  e 
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-1  K„ +  far  (?'ucT> 
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where  a  -  50/CrR_,  a„  =  0.26, 

1  t  V  l 


The  longitudinal  cross-correlation  obtained  from  Eq,  (9)  is 
plotted  in  Figure  (12) .  The  positive  isocorrelation  contours  are  symmetric 
about  the  line  C/U^j  *  constant.  Negative  isocorrelation  contours  are  also 
shown.  Although  the  negative  correlation  coefficients  are  small  compared 
to  the  positive,  the  negative  correlation  is  extended  over  a  larger 
separation  distance. 

For  a  maximum  correlation  with  an  optimum  time  delay  there  exists 
a  velocity  Uc  with  which  the  fluctuation  patterns  are  convected  down¬ 
stream.  For  large  separation  the  convection  velocity  is  obtained  from 
Eq.  (9)  by: 


3R(€,t) 


In  this  region  the  maximum-maximorum,  the  envelope  exp-|c|/a^6  is  quasi 
tangent  to  the  peak  at  ^max/T*  ^or  an  incompressible  flow  the  wall 
pressure  correlation  is  expected  to  closely  satisfy  the  boundary  conditions 
derived  by  Phillips  (1955)  and  Kraichnan  (1956),  that  is,  the  integral 


28 


ff  R(5,n,0)d€dn 

be  zero.  The  present  model,  however,  takes  into  account  the  compressi¬ 
bility  effect,  and  the  correlation  integral  does  not  vanish. 

From  the  spatial  and  temporal  extension  of  the  wall  pressure 
field  one  finds  it  useful  to  compute  the  integral  scale.  The  integrals 
are  defined: 


(11) 


Ap(0  =  2/  R(C)d5 

A  (n)  =  2/ R(n)dn 
p  ■&« 

and  A  (t)  =  2  f  R(t)dt 

p  t> 

Computation  and  comparison  is  given  in  the  table: 


(12) 
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3.03 

2100 
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0.11 

0.56 
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0.13 

Present 
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U  * 
e 

M  * 

£ 

1724,  U  = 
c 

1.98 

1356,  6  * 

.076 

0.52 

0.086 

0.11 

0.16 

The  longitudinal  scale  A(0/6  is  approximately  on  the  order  of 
the  boundary  layer  thickness  except  for  Bull's  result  which  is  one  order 


♦correction  applied  for  revision  of  data  by  author. 


Fig.  12.  Longitudinal  space  time  correlation  of  the  wall  pressure 
fluctuation  from  the  model. 
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of  magnitude  smaller.  The  lateral  scale  is  consistently  larger  than  the 

longitudinal  for  all  the  four  experimental  results.  The  ratio 

A  (5) /A  (n)  varies  from  one  tenth  for  supersonic  flows  to  about  one 
P  P 

sixth  tenth  for  the  subsonic.  It  is  of  interest  to  point  out  that  the  time 

integral  scale  and  the  ratio  of  the  spatial  scales  are  about  equal 

A  (5) /A  (n)  ss  U  A  (t)/6.  This  is  evidence  that  Taylor's  hypothesis  is 
P  P  e  p 

quite  valid  even  though  error  will  occur  in  the  extrapolation  of  the 
power  spectrum  it  (to)  to  zero  frequency. 

The  broad  band  spatial  cc  jlation  Fig.  (13)  shows  that  for  small 
spacing  the  correlation  field  is  nearly  isotropic  but  as  the  separation 
increases  the  longitudinal  correlation  becomes  negative  due  to  strong 
convection  while  the  lateral  decays  exponentially.  Figs.  (14  and  15)  show 
the  comparison  with  the  measurements  and  with  the  model  plotted  in 
terms  of  the  normalized  parameter  Ue5/<5Uc  (Eq.  7)  using  Taylor's 
hypothesis.  Agreement  is  quite  good. 

The  envelope  of  the  maximum-maximorum  is  plotted  in  Fig.  (16) 
using  the  equivalent  incompressible  Reynolds  number  C^Rg  to  compare 
correlation  made  at  different  Mach  numbers.  The  pattern  of  the  wall  pres¬ 
sure  confirms  that  the  increases  in  separation  reduces  the  higher  frequency 
contribution,  the  envelope  is  convex  for  larger  separation.  The  experimental 
results  show  some  scatter,  however  the  Reynolds  number  dependence  is 
consistent.  The  fact  that  the  convection  velocity  varies  with  separation 
distance  is  a  well-known  observation  which  Favre,  Gaviglio  and  Dumas 
(1967)  termed  celerity.  Typical  narrow  band  convection  velocity  is 
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given  in  Fig.  (17).  The  ratio  Uc(w)/Ue  almost  linearly  decreases  for 
large  io<S/Ue  ratio  as  indicated  in  some  previous  experiments. 

It  is  sometimes  convenient  to  work  in  the  frequency  domain 
rather  than  in  the  time  domain.  From  the  cross  correlation  Eq.  (9) 
one  obtains  the  cross  spectral  density 

-|cj/°t  6  -|n|/a,6  -itoS/U 

P(C,n;iw)  =  P(oj)  e  1  e  e  (13) 

0  ^  tu  5  +“ 


and  the  wave  number  spectrum 


P(KlfK2;u)  =  P(u>) 


l/ctjfi 


l/a2<s 


L^v)2  +  (ir  +  Ki)2JL(i/a26)2  +  (v 

as  03  tends  to  zero  the  wave  number  spectrum  varies  like  U.  Eq.  (14) 
indicates  that  the  maximum  occurs  when  9P(K^,u.<)/9K^  *  0,  correspondingly 
-w/K. 


(14) 


U 


1  max’ 

Fig.  (18)  shows  the  longitudinal  wave  number  spectra  obtained  from 
Eq.  (14).  The  contours  show  the  familiar  closed  loops  previously  reported 
by  Bradshaw  (1965)  in  a  low  speed  boundary  layer  flow.  The  main  difference 
being  that  the  flow  is  considered  to  be  semi-frozen  while  experimentally  one  can 
clearly  see  a  variation  in  convection  velocity  with  spatial  separation  and 
time  delay. 
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c)  Panel  Response 

The  panel  is  excited  by  a  steady  conveeted  supersonic  turbulent 
boundary  layer.  The  panel  deflection  is  small  compared  to  the  skin 
thickness.  It  is  assumed  that  the  response  is  linear,  stationary  and 
ergodic,  therefore  the  power  spectral  technique  is  employed  to  compute 
the  response,  i.e.,  Crocker,  Lyon  (1957),  Jacobs  and  Lagerquist  (1967), 
Maestrello  (1968),  Wilby  (1967)  and  White,  Cottis  (1968).  It  is 
assumed  that  the  plate  motion  does  not  modify  the  turbulent  boundary 
layer.  However,  the  forcing  function  should  contain  both  the  pressure 
over  a  rigid  wall  and  that  induced  by  the  panel's  own  motion.  In 
supersonic  flow  the  pressure  induced  by  the  panel  motion  cannot  be 
ignored  as  in  the  case  of  low  speed  flow,  because  the  supersonic 
convection  velocity  excites  strong  coupling  modes  to  the  acoustic  field. 

The  spectrum  response  is  mainly  controlled  by  acoustic  damping.  The 
pressure  due  to  the  panel  motion  is  comprised  of  waves  of  subsonic 
and  supersonic  velocities.  In  the  subsonic  region  the  damping  becomes 
purely  imaginary,  the  effect  is  to  increase  the  effective  mass  of  the 
panel  which  tends  to  reduce  the  resonant  frequencies  as  well  as  the 
response,  while  in  the  supersonic  region  the  damping  is  associated  with 
the  energy  radiating  field. 

The  equation  that  describes  the  pressure  on  the  surface  is  the 
well-known  Kirchhoff  solution.  Neglecting  the  viscous  term  and  the  terms 
that  are  nonlinear  in  the  surface  response  velocity  and  taking  into  account 
that  surface  pressure  is  twice  that  of  the  pressure  away  from  the  surface, 
and  applying  separately  the  real  and  image  of  the  boundary  layer  system 
from  Ffowcs  Williams  (1966),  one  can  write: 

P(x,t)  =  2T+  +  2P+  -  2P 

(15) 

0  =  2T_  -  2P+  -  2Py+ 


i&Wfm 
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where 
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Summing  the  real  and  image  pressure: 

P(x,t)  -  2(T+  +  T_)  -  4Pv+  (16) 

when  the  surface  Is  rigid 

Prigid(.i>t)  ‘  2(T+  +  T-)  (17> 

At  low  speed  flow  the  term  involving  the  pressure  due  to  the  panel 
motion  is  usually  neglected  since  the  mass  is  negligible  and  the  acoustic 
damping  is  small.  In  supersonic  flow,  the  panel  radiates  sound  in  the 
form  of  Mach  waves,  weaker  in  strength  due  to  inhomogeneity 
in  the  wave  number  components  but  with  loss  factors  sometimes  exceeding 
unity,  Maldanik  (1966).  The  measurements  made  of  the  displacement  spectra 
show  that  the  higher  frequency  modes  are  heavily  damped  unlike  the  panel 
excited  by  subsonic  flow.  Figure  (19)  shows  the  total  damping  factor  of 
a  12  x  6  x  .062  inch  titanium  panel  with  clamped  boundaries  excited  by  a 
turbulent  boundary  layer  at  Me  *  3.03.  The  bandwidth  Af  is  evaluated 
at  the  modal  half  power  point.  The  dotted  line  is  typical  of  experimental 
damping  for  a  panel  excited  by  low  speed  flow.  The  deviation  from  this 
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Fig.  19.  Damping  factor. 
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curve  at  higher  frequencies  is  attributed  to  the  wave  number  components 
becoming  supersonic.  The  effect  of  damping  tape  on  the  panel  is  not  as 
significant  as  might  be  expected  from  the  low  speed  case. 

Typical  displacement  spectral  densities  are  shown  in  Figs.  (20,  21  and  22) 

for  a  static  pressure  differential  of  0.06  psi  and  Fig.  (23)  for  a  static 

differential  of  14  psi.  For  the  spectra  taken  at  0.06  psi  the  panel 

is  enclosed  in  a  cavity  Fig.  (3)  while  the  spectrum  taken  at  14  psi  the 

panel  radiates  into  open  space.  As  expected  the  damping  factor  is  higher 

when  the  panel  radiates  into  open  space,  Strawderman  (1967),  Pretlove  (1964). 

The  broad-band  displacement  space  time  correlation  measured  along  the  panel 

centerline  in  the  direction  of  flow  for  pressure  differentials  of  .06 

and  14  psi  are  shown  in  Figs.  (24  and  25).  For  the  pressurized  panel 

the  response  is  mainly  due  to  coincidence.  The  ratio  between  the  wave  velocity 

U 

in  the  panel  and  that  of  the  turbulent  convection  velocity  is  =0.73 

c 

while  for  the  0.06  psi  pressure  differential  the  broad-band  cross 

correlation  shows  that  the  wave  velocity  of  the  panel  is  far  below  that 

U 

of  the  turbulent  convection  velocity  of  the  pressure  field,  =  0.29. 

c 

This  suggests  that  the  fundamental  mode  predominates  in  the  latter  case, 
contribution  also  results  from  the  modal  frequencies  of  the  unpressurized 
case  being  lower. 

The  contribution  of  the  closed  cavity  on  the  motion  of  the  panel  is 
small  when  the  depth  of  the  cavity  is  larger  than  one-half  the  larger  panel 

dimension,  Pretlove  (1964),  Dowell-Voss  (1963),  as  in  the  present  test. 

The  main  effect  is  on  the  damping  rather  than  the  panel  resonant  frequency, 
negligible  interaction  is  expected  between  the  fluid  and  the  plate.  In 
the  case  of  small  cavity  depths,  however,  the  principal  effects  on  the 
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Fig.  24.  Broad  band  space  time  correlation  of  panel  displacement. 
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fundamental  and  symmetrical  modes  will  be  that  of  an  aerodynamic  spring, 
whereas  for  the  antisymmetric  modes  only  the  aerodynamic  virtual  mass 
effects  exist. 

Fig.  (20)  shows  some  alteration  in  mode  shape,  for  example  the 
experimental  m  *  5  n  *  1  mode  is  strongly  coupled  to  the  m  =  3  n  =  2 
mode.  This  last  mode  should  be  zero  for  an  orthogonal  system  at  the 
measured  position.  A  similar  effect  is  seen  in  that  the  m  =  2  n  *  1 
and  m  =  4  n  =  1  are  also  non-zero. 

Fig.  (26)  shows  the  data  of  Fig.  (25)  replotted  on  the  £  -  U  t 

CP 

plane.  The  isocorrelation  contours  are  inclined  at  approximately  45  degrees 
indicating  a  convected  nature  and  the  decay  feature  of  the  correlation 
coefficient. 


AT  £  =  0  X=X'=  0.166ft. 
AT  Tj  —  0  Y  =  V'=  0.25ft. 

Ue=  2,100  ft. /sec. 


IV.  Analysis  of  the  Panel  Response 


The  response  of  linear  system  to  random  forcing  function  is 
obtained  utilizing  Green’s  functions.  The  system  (panel)  is  excited  by 
random  pressure  such  as  the  turbulent  boundary  layer,  therefore  the  assembled 
average  or  correlation  of  the  surface  displacement  is  expressed  in  terms  of 
the  correlation  of  the  pressure  field, 


_ _  .t  -t'  a  b  a  b 

y(x,y,t)y*(x',y’,t')  =  /  dt  dt 'Q  [  dx  f  dy  f  dxM  dy’ 

-03  -oo  J0  J0  J0 

g(x,y,t|Xo,yo,to)g*(x’,y’,t’|x;,y;,t;)<f(Xo,yo,t0)f*(x;,y;,t;)> 


(18) 


where  the  forcing  function  for  the  boundary  layer  pressure  field  is  given 
by  Equation  (9). 
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It  can  be  shown  (Dyer  1958)  for  a  plate  with  arbitrary  boundary  conditions, 
the  Green's  function  satisfies  the  classical  plate  equation  and  can  be 
represented  by: 
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where  $__(x»y)  are  the  normalized  eigenfunctions,  a  is  the  modal 
nm  inn 

damping  and  w  is  the  damped  resonance  frequency.  The  normal  mode  satis¬ 
fied  the  following  equation 


3Y 


B(l-in)V4Y  +  +  B  aT 

mn  .2  oT 


mn 


0 


(2 


3t 


where  B  is  the  bending  stiffness,  M  is  the  mass  per  unit  area  of  the 
plate,  B  is  the  viscous  damping  coefficient  and  n  the  .loss  factor 
associated  with  hysteretic  damping. 

The  modal  damping 
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where  T  is  the  eigenvalue  of  Eq.  (21). 

Equation  (13) involves  the  product  of  two  doubly  infinite  sums,  the 


cross-product  of  the  two  modes  (m,n)  and  (p,q)  is: 


f(-r-x)g.,(x)  +  e  f(T+x)g.(x) 


The  cross  spectral  density  is  obtained 
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In  computing  the  spectral  density  one  has  to  specify  the  eigenfunction 

a*0  % 

<}>(x,y)  and  <t>( - ,  - ),  Only  for  a  simply  supported  boundary  the  exact 

mi  nit 

solution  is  known,  and  opportunity  has  been  taken  of  this  solution  to  carry 
out  the  closed  form  solution  for  the  response  cross  spectral  density 
(appendix  b).  The  simply  supported  boundaries  are  not  physically  realistic; 
one  has  to  consider  spring  type  edges  conditioned  to  allow  rotation  and 
translation  constrained. 

Yang  (1950)  used  the  Ritz  method  to  find  the  dynamic  characteristics 
of  rectangular  plates  with  continuum  spring  boundary  conditions.  More 
recently  this  method  has  been  brought  up  by  Izzo,  Baylor,  and  Robideaus 
(1960) ,  The  present  analysis  does  not  restrict  the  usage  of  any  particular 
method  (i.e.,  Ritz,  Galerkin  or  Rayleigh)  for  example  in  the  case  of  the 
present  experiment  the  panel  boundaries  are  clamped;  Warbur ton’s  (1954) 
approximation  is  used  in  the  computation. 

The  acoustic  power  radiated  by  the  panel  is  obtained  by  the  same 
procedure  used  to  compute  the  acoustic  power  radiated  by  the  turbulent 
boundary  layer  described  in  the  next  section.  Using  this  result  expressed  in 
terms  of  the  panel  displacement  cross  spectral  density  Y(X;w)  rather  than 
the  pressure  cross  spectral  density  over  the  rigid  wall  P(X;u)  the 
radiated  power  spectral  density  is  given  by 
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where  A  »  /(x-x')^  -  (y-y’)^  . 
~  '  o'  w 
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The  pressure  radiating  modes  of  the  flexible  panel  mounted  on  an  infinite 
rigid  plane  is  in  phase  with  panel  velocity,  the  wave  motion  has  supersonic 
phase  velocity,  the  energy  losses  appear  as  a  damping  of  the  panel  response. 
While  the  subsonic  components  (near  the  plate  surface)  appear  as  additional 
mass,  they  are  of  predominant  lower  frequency. 
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V.  Comparison  with  Measurements 

A  good  comparison  is  made  with  calculation  for  two  panel  configurations, 
one  with  0.06  psi  pressure  differential  enclosed  in  a  cavity  (Figs.  20,  21 
and  22)  and  the  other  with  14  psi  pressure  differential  radiating  into  open 
air  (Fig.  23).  The  calculations  are  based  on  the  analysis  discussed  in  the 
previous  section  which  does  not  account  for  the  additional  forcing  field,  the 
pressure  induced  by  panel  motion.  However,  the  damping  used  in  the  calculation 
are  those  values  obtained  experimentally.  As  indicated  previously,  the  damping 
is  mainly  acoustic  and  thus  one  can  ignore  the  mechanical  damping.  The  subsonic 
component  of  the  velocity  is  out  of  phase  with  panel  velocity  and  therefore  no 
dissipation  of  energy  occurs.  The  rise  of  the  virtual  mass  due  to  the  subsonic 
component  is  small  for  an  air  medium  and  it  can  be  ignored.  By  making 
e  *  1  in  the  Green's  function  and  taking  into  account  the  pressure 

induced  by  the  panel  motion  the  acoustic  damping  will  be  a  part  of  the  solution. 
This  calculation  procedure  is  a  result  of  the  present  work  and  will  soon  be 
reported.  It  is  expected  that  the  new  method  (Maestrello,  McDaid)  in  computing 
the  response  may  not  differ  more  than  a  factor  of  two  (in  the  present  case). 

In  addition  to  the  case  where  the  plate  motion  alters  the  turbulent  boundary 
layer  a  perturbation  pressure  similar  to  Kuo  (1965)  in  study  of  the  surface 
roughness  has  to  be  added. 

The  pressurized  panel  has  lower  displacement  spectra  than  the 

2  2 
3  Y  a  Y 

unpressurized  one.  In  fact  the  static  loading  N  — r  +  N  — r  +  AP  is 

X  3x  y  3y 

the  contribution  of  it.  The  shift  in  frequency  from  the  pressure  of 
0.06  to  14  psi  is  from  475  to  595  in  the  fundamental  mode.  This 
is  accomplished  by  a  change  in  static  deflection  of  about  0.07  inch  at 


the  panel  center,  effectively  Increasing  the  stiffness  of  the  panel. 
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structure  by  means  of  a  statistical  method.  In  a  similar  fashion  utilizing 
the  temporal  and  spatial  correlation  of  the  wall  pressure,  separation  is 
obtained  of  the  acoustic  mode,  that  is  the  weakest  mode  within  the  flow 
field,  from  the  vorticity  and  entropy  modes  being  the  stronger  by  far. 

The  acoustic  power  predicted  is  obtained  by  imaging  techniques  where 
the  properties  of  the  far  field  noise  are  connected  to  the  statistical 
properties  of  the  boundary  layer  via  near  field  pressure  over  the  wall. 

Powell  (i960)  showed  that  a  plane  rigid  surface  reflects  the  sound  radiated 
by  an  equivalent  quaarupole  located  in  the  flow  stream  above  the  surface. 

The  surface  imaging  pressure  field  is  identical  to  measurements  made  of 
space-time  correlation  of  the  double  time  derivative  of  the  Reynolds  stresses 
ai.d  integrated  over  the  turbulent  field  (Lighthill  results).  Where  the  flow 


4  I 

lif 


is  over  the  surface  the  reflection  from  the  surface  could  be  distorted  due  to 
the  action  of  viscous  stress;  it  has  been  argued,  however,  that  the  normal 
stress  on  the  surface  is  negligible*.  With  this  assumption  the  problem  is 


equivalent  to  that  of  an  image  boundary  layer  in  which  the  acoustic  power  radiated 
per  unit  volume  of  turbulence  could  be  related  to  the  correlation  area  of  the 
pressure  on  the  surface.  Since  the  surface  is  rigid  no  acoustic  energy  can  be 
generated  but  it  is  generated  within  the  turbulence  due  to  variation  of  the  eddy 
strength  caused  by  the  fluctuating  volume. 

Laufer  (1961)  made  the  first  measurements  of  sound  radiated  by  a 
supersonic  turbulent  boundary  layer,  by  means  of  a  hot  wire  annemometer, 

*Leehey  (1968)  recognized  that  the  neglected  surface  integral  in  the  evaluation 
of  the  wall  pressure  fluctuation  gives  rise  to  oscillatory  wall  shear  stresses. 
Leehey's  observation  came  to  the  author's  attention  after  this  manuscript  was 
written,  therefore  the  results  of  the  implication  above  are  not  yet  known. 
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correlating  in  space  and  time  outside  the  boundary  layer.  The  loss  of 
energy  from  the  turbulence  shearing  stress  by  acoustic  radiation  is  very 
small,  therefore  the  detection  and  measurement  made  by  Laufer  was  quite  a 
challenge.  His  results  agree  qualitatively  with  the  model  of  the  radiation 
field  as  predicted  by  Phillips  at  higher  Mach  numbers. 

Wilson  (1960)  also  measured  the  sound  radiated  from  a  subsonic 
boundary  layer  around  a  rotating  cylinder.  He  placed  the  rotating  cylinder 
into  a  reverberation  room;  the  acoustic  pressure  radiated  from  the  boundary 
layer  was  obtained  after  eliminating  the  mechanical  noise  of  the  rotating 
mechanism. 

Ffowcs  Williams  and  Maidanik  (1965)  computed  the  mean  square  radiated 
pressure  in  terms  of  surface  pressure  fluctuations  using  Lighthill 
formalization.  Their  results  well  agree  with  the  measurements  made  by 
Laufer. 

In  formulating  the  present  problem,  let  us  recall  some  fundamental 
flow  field  equations.  Assuming  the  flow  is  incompressible,  uniform 
density  p,  the  pressure  p  obeys  the  following  relationship 

V2  P(x)  =  -  f(x)  (29) 

the  free  field  solution  f  *  0  at  infinity  is  given  by 


or 


(30) 


(31) 


the  surface  being  infinite  rigid  and  plane  an  additional  boundary  condition 
is  imposed. 
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Then 
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from  Gauss* s  theorem 
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The  general  solution  is  obtained  by  adding  Equations  (31  and  33) 

p(ii=^/^f  *  “  £  /[t^  If  - p  s  fer)]  * 

is  the  solution  of  Poisson's  equation  for  the  case  of  a  rigid  boundary. 

If  one  takes  the  pressure  P  at  the  retarded  time  |x-y|/c  where  c  is  the 

propagation  velocity,  the  above  equation  becomes  the  solution  of  the  in- 
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homogeneous  wave  equation.  Evaluating  —  one  obtains: 
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Substituting  into  (34) 
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The  above  equation  is  the  well-known  Kirchhoff  solution.  The  volume 
integral  expresses  the  contribution  of  sound  from  the  turbulence,  while 
the  surface  integral  is  the  contribution  of  the  sound  upon  the  surface. 
By  imaging,  the  two  integrals  are  identical.  Three  terms  appear  in 
the  surface  integral,  the  first  gives  the  contribution  of  the  usual 
pressure  on  the  surface,  the  second  describes  the  radiated  pressure 
from  the  apparent  sources,  and  the  third  term  describes  the  near  field 
pressure  which  decays  1/r  faster  than  the  radiated  pressure 
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The  corresponding  pressure  correlation 
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From  the  theory  of  Fourier  transforms  one  obtains 


R(x,y;u)=  -ai--  '2°2  f  f 
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The  power  spectral  density  per  unit  area  of  the  acoustic  pressure 
field  is  obtained  by 
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where  P(A^u)  is  the  power  spectral  density  of  the  pressure  field. 

The  integrand  behavior  is  such  that  as  u|l|/c  tends  to  zero  the  term 
{  }  — ~  ss  u  /c  ,  meaning  the  radiated  spectrum  will  peak  higher 

IM 

than  the  near  field  pressure,  a  result  consistent  with  Lyamshev  (1961). 
For  turbulent  boundary  layer  a  model  of  P(A;u)  is  given  by 
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Fig.  27.  Acoustic  power  level  from  turbulent  boundary  layer  over  a 
rotating  cylinder. 


Wilson's  experiment  is  used  for  comparison  with  the  present  analysis 
Fig.  (27).  In  order  to  compute  the  acoustic  power  radiated  by  the 
turbulent  boundary  layer  on  a  rotating  cylinder  one  needs  to  evaluate 
the  cross  spectral  density  P(^,aj).  The  model  of  the  turbulent 
boundary  over  a  flat  plate  is  adopted  for  PCX^w)  and  Wilson's  near 
tield  pressure  spectrum  is  used  since  at  lower  frequencies  the  near 
field  power  spectrum  differs  from  that  over  a  flat  plate.  The  result 
ing  comparison  is  satisfactory. 

For  subsonic  flow  the  cross  spectral  density  P(X,m)  is 


proportional  to 
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eK  e 


and  the  acoustic  power  is  proportional  to 
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while  for  Mach  wave  radiation  the  cross  spectral  density  P(X,u>)  is 


proportional  to 


2  3 
p  U  6 
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and  the  acoustic  power  is  proportional  to 
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In  attempting  to  compute  Laufer's  Mach  wave  radiation,  a 


difficulty  is  encountered  since  he  measured  the  intensity  at  a  point 


r 


while  this  analysis  obtains  the  power  per  unit  area.  In  the  comparison 

/M  \ 

one  would  have  to  account  for  the  Doppler  shift  (  —•  -  1 J  since  Laufer's 
measurements  are  in  a  moving  stream  along  the  Mach  line.  The  condition 
for  Mach  wave  radiation  is  such  that  |Mc  -  M|  >1.  This  implies  that 
the  comparison  has  to  be  made  for  Mg  >  3.5. 

Phillips  has  shown  that  radiated  sound  arises  from  eddy  Mach 
waves  which  are  generated  by  some  wave  numbers  of  the  turbulence  in 
those  layers  for  which  the  difference  between  the  mean  velocity  of  the 
fluid  outside  and  the  local  eddy  velocity  is  greater  than  the  speed  of 
sound.  This  mechanism  appears  consistent  with  the  present  results. 


APPENDIX  A 


Prove  that  the  displacement  cross  correlation  is  even  in  t 
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Let  us  define 
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APPENDIX  B 

The  following  is  the  closed  form  solution  of  the  panel  deflection 
cross  power  spectral  density  (Eq.  24)  for  simply  supported  boundaries 
without  cross  modal  coupling.  From  Maestrello  (1967)  one  can  obtain  the 
following  cross  correlation  of  the  panel  displacement. 
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I  is  an  even  function  of  t  and  Y(x,y,t)  Y  (x’jy’jt’) 
is  also  even  as  shown  in  (appendix  a);  one  can  obtain  the 
auto  correlation  by  setting  x  *  x'  and  y  *  y'  in  (IB). 
The  displacement  cross  spectral  density  is  obtained: 
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From  (2B)  and  integrating  with  respect  to  x  one  obtains: 
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Simplify  the  above  equation  by  integrating  out  the  g(x)  and 
f(z)  integrals; 
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Finally  the  f(y)  integral  reduces  to: 
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The  cross  spectral  density  is  finally  reduced  to: 
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